# -*- coding: utf-8 -*-
from scipy import *
from plate import Plate
from draw_func import *
import time

if __name__ == '__main__':
	
	beta, b = 1.,100.
	a = b*beta
	gamma = 0.0
	M,N = 10, 10
	
	p = Plate()
	p.M,p.N = M,N
	
	Beta = linspace(0.4,3.0,100)
	Gamma = [0.0,1.0,10.0,50.0]
	for gamma in Gamma:
		p.gamma = gamma
		K1 = []
		for beta in Beta:
			p.beta = beta
			p.bucklingw()
			K1.append(abs(p.K[0]))
		
		pfigure(Beta,K1,
			mk = '',
			pL='$ \\gamma_{sf} = %.1f $'%gamma)
	#~ show()
	#~ p.gamma = gamma
	#~ p.BC = 2
	#~ p.ds = 0
	
	#~ t0 = time.clock()
	#~ p.buckling()
	#~ t1 = time.clock()
	#~ K,A = p.K,p.A
	
	
